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bX), We analyze quantum mechanical systems using the non-perturbative renormalization 

group (NPRG). The NPRG method enables us to calculate quantum corrections systemati- 
■ cally and is very effective for studying non-perturbative dynamics. We start with anharmonic 

oscillators and proceed to asymmetric double well potentials, supersymmetric quantum me- 
OO , chanics and many particle systems. 

(N . 

t-h ' §1. Introduction 

> : 

f- — , , The non-perturbative renormalization group (NPRG) has been formulated through 

t-h ' analyses of critical phenomenal and applied to non-perturbative studies of statis- 

tical mechanics and quantum field theories. It has been established as a powerful 
, tool for analyses of non-perturbative dynamics in systems of many (infinite) degrees 

of freedom, because it allows for the evaluation of fluctuations without recourse to 
perturbation series. Several types of non-perturbative {exact) renormalization group 
Q-i, equations have been derived by integration with scale decomposition and have been 

■y \ applied to various systems. In this article, we apply the NPRG 

method to quantum mechanical systems, that is, systems of finitely many degrees of 
^ freedom B'Q^'tO) to analyze their non-perturbative dynamics. 

Generally, there are two types of non-perturbative quantities. One corresponds 
to the summation ojLall orders of a perturbative series, which could be related to 
Borel resummation.lHP The other is an essential singularity with respect to a coupling 

^ i constant A , which has a structure like e"^.E§'H We are not able to expand such 

a singular contribution around Ao = 0. A singularity of this type appears in the case 

of quantum tunneling. For example, in a symmetric double well system, there are 

two degenerate energy levels at each minimum, which are mixed through tunneling 

i_ 

to generate an energy gap AE ~ e A o . The exponential factor is known to result 
from the free energy of topological configurations, i.e., instantons. 

In this article, we first summarize how to analyze quantum mechanical systems 
using the concept of NPRG and check to what extent NPRG can be used to evaluate 
non-perturbative effects quantitatively. The NPRG equation we employ dusre is a 
local potential approximated Wegner- Houghton (LPA W-H) equation, £3>>la> which 
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we use to analyze quantum anharmonic oscillators and asymmetric double well sys- 
tems. In contrast to the symmetric double well system, the standard instanton 
method does not work for an asymmetric potential, and the much more sophisti- 
cated method of the valley instanton has been developed for their treatment. H'ElP 
The NPRG method is found to work for asymmetric potentials as well as for sym- 
metric potentials, because NPRG does not rely on parity symmetry. 

We proceed to analyses of more complicated systems, super symmetric quantum 
mechanics (SUSY QM) and many particle systems. SUSY QM is a toy model for 
dynamical SUSY breaking. Ef'E^ Although, in general, there is no spontaneous sym- 
metry breaking in systems with finitely many degrees of freedom, some extraordinary 
symmetries, such as SUSY, can be broken even in quantum mechanics. SUSY break- 
ing is a highly non-perturbative phenomenon because of the non-renormalization the- 
orem, and we will see that NPRG should be applicable for non-perturbative SUSY 
breaking. 

In addition, analyses of quantum many particle systems have become very im- 
portant with recent developments in nano-technology. Solving the problem of how 
the quantum coherence of a variable of a target system is affected by other variables 
(the environment) is quite important. For example, it is necessary for realization of 
qubit for quantum computers. However, standard methods that are well suited for 
treating systems of one degree of freedom, the Schrodinger equation, instanton, etc., 
do not work well in such complicated systems. We believe that NPRG is versatile 
enough to analyze such systems. As a first step, we analyze quantum tunneling 
phenomena in two particle quantum systems. 



In this section, we briefly summarize the formulation of NPRG with D-dimensional 
real scalar field theory. 

2.1. Scale decomposition 

In the NPRG method, the theory is defined by the Wilsonian effective action 
<S/i [0]. This is an effective theory with an ultraviolet energy cutoff A: 



We decompose the path integration variable 4>(p) into two parts with respect to the 
momentum scale p as 



§2. Non-perturbative renormalization group 




(1) 




< \p\ < A — AA : lower modes, 
A - AA < \p\ < A : shell modes, 



and transform the partition function Z as follows: 
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V(f) < e 
V<j> 



-Sa[4>< 



— ■? mt \ 



[</><] g-AS.4 [</><] 



-5a-z1/i[</'<] 



where 



Z\Sa[0<] = -log 



c in 

e ^ 



= / V<t> a 



-Sa[4>s 



e b A 



(2) 



(3) 



We understand the shell mode path integral (• • •} ^ as the renormalization trans- 
formation. If we evaluate it by perturbative expansion with respect to coupling 



constants, we obtain the so-called perturbative renormalization group equations. liSP 
Of course such equations are valid only in the weak coupling limit. Instead, we take 
the limit AA — > to define NPRG equation, which is the fundamental procedure. 

2.2. Derivation of the NPRG equation 

Taking the limit AA — ► 0, we can express the renormalization transformation as 
a differential equation, 

9Sa ,. Sa — Sa-aa 1 



lim 

dA AA-*o AA 
Then, we can expand Sa[<P] as power series in 

SS A 



lim — — - log 

AA-+0 AA 



_ c in 

e a A 



(4) 



S A [<t>] = S A 



Js] 



1 

+ 2 



shell 



shell 0~4>{P) 
(P) 



(P) 



S =O 



S<f>(p)6</>(q) 



,(q) + 0(AA 2 



(5) 



</> s =o 



Since an O(AA) calculation is required for the evaluation of the derivative (4), the 
shell mode path integral (• • can be evaluated exactly using a Gaussian integra- 
tion. Then, the fundamental differential equation is derived as 



A 



dS A 
dA 



shell 



log 



5 2 S A 



+ 



SS A 



S 2 S A 



SS a 



• (6) 



This is known as the Wegner-Houghton equation. @) It represents exactly the cutoff 
A dependence of the Wilsonian effective action Sa- Its right-hand side is generally 
called a (3 functional. 

2.3. Approximations 

Although the Wegner-Houghton equation is exact, we cannot solve it without 
some approximation in practice. In this article we employ the local potential approx- 
imation (LPA), which means that we ignore corrections to derivative interactions. 
It can be considered the leading order of the derivative expansion of Sa, 



Sa[<P] 



d D x 



VA[<t>} + ~K A [<t>]d^ 



+ 



(7) 
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To make this approximation, we substitute the zero mode for the lower mode (p < (p): 

Mp)^v(^) D 5 d ( P ). (8) 

The local potential approximated Wegner-Houghton (LPA W-H) equation is then 
obtained as follows: 



A 



9Va 
dA 

A D 



-^A D log (l + 
Jdf2 D 



1 d 2 V A 
A 2 dp 2 



0) 
(10) 



(2vr) D " 

This is a two-dimensional partial differential equation for Va (p)- 0^ Its right-hand 
side is called a (3 function. We solve it mainly using numerical methods. 

To obtain an intuitive understanding of this equation, we proceed to further 
approximation, the operator expansion. We expand V e g as power series in ip: 



N 



a n (A) 



n=0 



-<p 



(11) 



The partial differential equation is then reduced to a set of ordinary differential 
equations for the coupling constants {a n (A)}: 



dao 
~dA 
dai 
~d~A 
da 2 
~dA 

da-s 
dA 

da^ 
~dA 



A 



«2 



A 2 



= -^A° 
2 


«3 


.A 2 + a 2 


= -^A° 
2 




A 2 + a 2 


= -^A D 
2 


a 5 


A 2 + a 2 


= -^A D 
2 


a 6 


A 2 + a 2 



+ 



2a| 



4a 5 a 3 



{A 2 + o 2 ) 



(12) 
(13) 
(14) 

(15) 



3a^ 



+ 



(A 2 + a 2 ) 
12a4a| 



6aq 



(A 2 + a 2 f (A 2 + a 2 ) 



(16) 



If the results of these equations converge as the order of the truncation, N, becomes 
large, we regard them as solutions of the LPA W-H equation. HIP This expansion 
allows us to treat differential equations more easily and to understand the origin 
and structure of quantum corrections physically. For example, a correction to odd 
n couplings cannot be generated from even n couplings only. This implies that if we 
choose an initial potential Vo(ip) as Z 2 ((p *-* —ip) symmetric, the solutions of the RG 
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equations flow in Z2 symmetric subspace; that is, the NPRG equation does not break 
the global symmetry of the system. However, we should note that in some extreme 
cases, the operator expansion is not good, and leads to pathological behavior, as 
seen in §4. 

The constant part of Va, ao, is given by the vacuum bubble diagramso and 
is usually ignored. However, we retain it here, because it plays a crucial role in 
supersymmetric theories. 



§3. NPRG analysis of quantum mechanical systems 



Making use of the LPA W-H equation, we analyze systems in quantum mechan- 
ics, which is D = 1 real scalar theory with a single dynamical variable, x(r).H2P 

3.1. Physical quantities 

The LPA W-H equation for quantum mechanics is given as follows: 

4&~5F'*( 1 + ™)- 

We solve it by lowering A from the initial cutoff Aq, where the initial potential Va 
is given by the potential term Vq(x) in the original action, 



S[x] = Jdr I Ix 2 + y (x)|. (18) 



In the infrared limit A — * 0, we obtain the effective potential V e ff(x) = lini/t^o Va(x) 
from which physical quantities are evaluated. 

First, the expectation value of x in the ground state \ f2), 

(x) = (n\x\n), (19) 

is determined by the stationarity condition, 

dV cS 



= 0. (20) 

x=(x) 



dx 

The ground state energy of the quantum system is given by 

E = {f2\H\f2) = V eS (x = {x)). (21) 

Also, we obtain the energy gap of the system through the following expressions of 
the two-point correlation function: 

<™*(o)i«> = i d i^t E , +i ":_ Ea) —^'-^^) 

<fl| T H r )m \0) ^ I f^Vj^ = 5^—*, (23) 



* For other types of NTpRG equations, such as the Legendre flow equation, which is derived by 
means of a cutoff function,^ we cannot evaluate the vacuum bubble diagrams properly without 
the prescription of subtracting the contribution of the cutoff function from the constant part of Va- 
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where D„ = 2\C n \ 2 (E n - E ), C n = (n\x(Q)\f2), J2n D n = 1 and m 2 eS is the 
curvature of the effective potential at the minimum. Comparing the damping factors 
of (22) and (23) in the r — > oo region, the energy gap AE = E\ — Eq is obtained as 
follows :□ 



AE = m eS 



d 2 V r 



eff 



dx 2 



(24) 

x=(x) 



Generally, the NPRG method does not evaluate the wave functions of the tar- 
get system directly. Instead, it yields information regarding Green functions. For 
example, the n-th moment of x corresponds to information concerning the ground 
state wave function ipo(x): 

M n = (f2\ x n \Q) = fdx x n \Mx)\ 2 ■ (25) 



We now give some examples of calculations of M n in a ^-symmetric (namely, (x) = 
0) system. The two-point function M<i is calculated as 

m 2 = { n\ x 2 \n )c LP = A I f^^r = (26) 

J 2tt E z + m^ ff 2m eS 

where the subscript c denotes a connected function. In a similar way, the four-point 
function M4, 



M 4 = (J?| x 4 \f2) c + 3M 2 2 , (27) 

|4 v cff 



is calculated by means of the LPA four-point coupling A e fj = '' 



,n\ -4m\ LPA a I \ f dEidE 2 dEz 
(f2\ x \n) c = -4! A eff J 



El+m 2 cS (E 2 -E x f+ml s {E 3 -E 2 ) 2 + m 2 cS E 2 +m\ 



2 ' 

eff 



(28) 



In this way, we are able to connect the effective couplings obtained from the NPRG 
with M n obtained from the ground state wave function. All of this information 
concerning the ground state taken together provides full information regarding the 
quantum system, including excited states. 

3.2. Example: harmonic oscillator 

To illustrate the characteristics of the NPRG analysis, we now consider the 
harmonic oscillator. The initial potential is chosen as Vo(x) = ^m 2 x 2 at the initial 
cutoff scale Aq. In this case, we can solve the LPA W-H equation analytically, and 



* For further details with regard to this relation, see Appendix A. 
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we obtain 



ao(A) = a (A ) + 



2vr 



plog 



1+p 2 
p 2 



+ 2 tan l p 



V a 2( A o) 



V /a 2( yl o) 



(29) 



(30) 



02(A) = a 2 (A ). 

Since the initial conditions are (an(A))> 02(A))) = (0, m 2 ), if we take the simultaneous 
limit Aq — > 00, A —> 0, (ao(A), 02(A) = (f ,m 2 ) is obtained. Although a2(A is free 
from quantum corrections and does not run, oq runs and produces a zero-point energy 



2 • 

Figure.[y plots the actual running of 
ao and shows that it is limited to a fi- 
nite energy region that depends on the 
mass scale, m. Ultraviolet finiteness is 
a typical feature of quantum mechanical 
systems, and it implies that the theory 
is finite, even in the Aq — > 00 limit. Con- 
trastingly, the infrared finiteness in Fig- 
ure.|l] is related to the decoupling prop- 
erty that a heavy particle cannot prop- 
agate in the low energy region. Such 
ultraviolet finiteness and infrared finite- 
ness enable us to obtain physical quanti- 
ties even through numerical calculation 
within a finite energy scale region. 



< 

c 
CO 



Running of Vacuum Energy 

harmonic oscillator 



5.0 — 



4.0 
3.0 
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Fig. 1. Running of ao 
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§4. Analysis of anharmonic oscillators and double well systems 

4.1. Symmetric single-well potential 

Now we proceed to analyze quantum mechanical anharmonic oscillators and 
double well systems. First, we consider a symmetric single-well potential, 

V (x)= X Q x 4 + l -x 2 . (31) 

Our interest is to compare our NPRG results with the perturbative series. First, the 
LPA W-H equation (17) is solved numerically, and we thereby obtain an effective 
potential V e s. The flow of Va is shown in Figure.|2[ Quantum corrections raise the 
potential and make its slope steeper. In Figure]^, we display the energy spectrum 
calculated with the relations (21) and (24). We refer to the results obtained by a 
numerical analysis of the Schrodinger equation as the "exact results." 
The perturbative series of E n is the asymptotic series 

En = [n + -] + jAopn 2 + 2n + 1] - ^A^[34n 3 + 51n 2 + 59n + 21] + • • • . (32) 
2 4 o 

It diverges even in the weak coupling region. Note that the Borel resummation of 

the perturbative series works well in this case, and gives quantitatively good values. 
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Change of Potential Energy Spectrum 

LPAW-H, VO-1 single well 




Fig. 4. M 2 = (n\ x 2 \Q). Fig. 5. Effective coupling A cff . 



However, even in the lowest order approximation (LPA), the W-H equation can 
evaluate the energy spectrum almost perfectly. Therefore, we conclude that the 
NPRG does sum up all orders of the perturbative series in the correct manner. 

We also display the two-point function M2 = {f2\ x 2 \£2) in Figure]!] and the 
effective four-point coupling constant A e g in Figure.^. In both of these cases as well, 
the NPRG results give an almost perfect fit. 

4.2. Symmetric double-well potential 

Next, we consider the ^-symmetric x «-> — x double-well potential, 

V (x) = A x 4 - X -x 2 . (33) 

In quantum mechanical systems, this Zi symmetry never breaks spontaneously, be- 
cause the x mode tunnels through the potential barrier, and the ground state is 
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uniquely realized. In fact, in the NPRG evolution of the effective potential, the ini- 
tial double-well potential finally becomes a single well, and an energy gap (effective 
mass) arises (Figure.^). 

In this system, there is no well-defined perturbation theory. A standard tech- 
nique to obtain the energy gap AE is the dilute gas instanton calculation. This is a 
semi-classical method based on the one-instanton solutions 



1 



x c \{T) 



: tanh 



The one-instanton contribution to the partition function Z is 

-1/2 



(34) 



Z ~ Te~ s[Xc 



S[x cl ] 
2vr 



det' ( 



X x c i 



det 



5 2 S 

d.r 2 



x=±- 



= T- 



(35) 



where T is an imaginary time volume. Assuming that instantons do not interact 
with each other, we can evaluate the multi-instanton contribution to Z (the dilute 
gas instanton approximation), and we obtain the energy gap 



AE = An 



h\[2 

7rAo 



(36) 



which has the structure of an essential singularity originating from the one-instanton 
action. The singularity coefficient obtained from the instanton method is known to 
be exact in the vanishing Aq limit. 1 



fL9h 



Change of Potential 

LPA W-H , V=0.1 



Energy Gap 

double well 



> 




-1.0 



-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 



Fig. 6. Flow of the potential. 




0.40 



Fig. 7. Energy gap estimates. 



In Figure.^ we display the energy gap evaluated using various methods. The 
NPRG results are very good in the strong coupling region, while the perturbation 
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cannot be applied in this double-well system, and the dilute gas instanton method 
is not at all effective, because it is valid only in the very weak coupling region.u 
Therefore, the NPRG method should provide a powerful tool for the analysis of 
tunneling, at least in such regions. However, our NPRG results deviate from the 
exact values as Ao — ► 0, which corresponds to a very deep well. Because the (3- 
function becomes singular in this region, the NPRG results become unreliable. We 
believe that the cause of the difficulty is the LPA approximation scheme that we 
adopt. It is important to note that the respective coupling regions in which the LPA 
W-H equation and the dilute gas instanton work well are_separated, and therefore 
these two methods should be regarded as complementary. t2l)'E3) 

We display the two-point function M2 = {Q\ x 2 \f2) in Figure]^ and the effective 
four-point coupling A e g in Figure.^. As in the case of the energy gap, the NPRG 
results are excellent here, except in the extremely weak coupling region. 

4.3. Flow diagrams 

We now more carefully consider the difficulty arising in the weak coupling region 
for the double-well potential. We employ the operator expansion (11) and investi- 

n + 2 

gate the flows of the dimensionless coupling constants a n = a n A — — . The flow 
diagrams elucidate the phase structure of the system. We display the flow diagrams 
for the N = 4,6,10 truncated potentials (Figure .|i~0|, Figure. [ll], Figure. 12) and for 
the potential without an operator expansion (Figure.|l~3|) . 

These flow diagrams reveal that the phase structure of the theory with a trun- 
cated potential (JV = 4, 6, 10) is somewhat strange. As mentioned above, there is no 



* It has long been known that the strong coupling expansion has a finite radius of convergence. 
Recently, variationaL-oectjurbation theory has become highly developed, and very accurate results 
have been obtained. Ej'Ex The region of coupling constant values in which these approaches are 
good is estimated as Ao > 0.08, which is almost coincident with the reliable region for our method. 
To elucidate the correspondence between the NPRG method and this improved perturbation theory 
is interesting. 



Non-Perturbative Renormalization 



N=4 



4.0 



F.P. /^f 


mil 




/ Broken 
/ Phase \ 





-1.0 0.0 



a 2 (dimensionless) 

Fig. 10. N = 4 truncation. 



N=10 




-1.0 0.0 



a 2 (dimensionless) 

Fig. 12. N = 10 truncation. 
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N=6 




-1.0 0.0 



a 2 (dimensionless) 

Fig. 11. TV = 6 truncation. 



No truncation 




a 2 (dimensionless) 



Fig. 13. No truncation. 



spontaneous symmetry breaking in these quantum mechanical systems. When we 
truncate the potential as in (11), there appears a non-actual fixed point and false bro- 
ken phases (Figures [l(| [II] and 12). The flow starting from the weak coupling region 
(Aq — > i.e. 04 — > 0) tends to be captured by the false broken phase, and we cannot 
obtain the correct result > 0. The region of the false broken phase becomes 
smaller as N becomes larger, and then for the LPA exact (no truncation) calculation, 
the false broken phase disappears (Figure]!^). However, even in the no truncation 
case, we cannot obtain reliable results for the flows that start from the weak coupling 
region, because singular behavior of the flow in the region near 0,2 = — 1 leads to large 
numerical errors. The results for the LPA W-H equation in Figure]?] were obtained 
from numerical integration of the partial differential equation without any trunca- 
tion.r§ 
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Energy Gap 

double well 



4.4. Other methods 



0.8 



0.6 



— exact result 
LPA W-H 
• LPA Legendre k=0 

LPA Proper Time RG 
i NLO Proper Time RG 
■ Density Matrix RG 



3 0.4 



LU 

< 



0.2 



0.0 



The NPRG equation we employ 
here is that with the local potential ap- 
proximation. The results in the weak 
coupling region can be improved by up- 
grading the approximation. The LPA is 
the lowest order of the derivative expan- 
sion, and a higher-order calculation can 
be carried out.E^ In this quantum me- 
chanical system, the second-order cal- 
culation of the Legendre flow equation 
does not improve the weak coupling re- 
sults. However, an analysis using the 
proper time renormalization group im- 
proves the LPA results considerably. c£> 
Also, although it differs from the NPRG 

methods in its formulation, the density matrix renormalization group is useful for 
this system. O We exhibit in Figure.[l4] tharesults for the energy gap obtained with 
various renormalization group approaches .□ 

4.5. Asymmetric double-well potential 

We proceed to consider the ^-asymmetric double-well potential 



0.0 



0.1 



0.2 



Fig. 14. Various RG results. 



Vq(x) = X x 4 - ^x 2 + hox, 



(37) 



where the linear term Hqx breaks the Z2 symmetry explicitly. In this system there 
are a stable minimum and an unstable minimum. 

How do we deal with the effect of such an asymmetric term? The NPRG method 
can treat this system in a manner that is quite similar to that for the symmetric 
system; it just changes the initial potential, while the LPA W-H equation does not 
change. Furthermore, when we apply the operator expansion (11), the situation 
becomes even simpler. The additional h$x term does not affect the running of other 
coupling constants, because the term 



dr hox(r) = hox(E = 0), 



(38) 



consists entirely of the zero energy mode, and generates no quantum corrections. 
Therefore, the NPRG equations for the coupling constants are the same as those in 
the symmetric case. 



* As for general non-perturbative methods, the auxiliary field method works very well both 
in the weak coupling and strong coupling regions. E3' Also, various improved perturbation theories 

S^e-beep .applied to the anharmonic oscillator and the double well system, giving similar results. 
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By contrast, the standard instanton method cannot be applied to such an asym- 



metric system, because the term 



S 2 S 
SxSx 



X — X c \ 



in (35) has a negative eigenvalue in this 



case. For actually unstable systems, this negative eigenvalue is converted to a decay 
rate for the system. This is a typical prescription for the 'bounce solution' calcula- 
tion. However, in the case of the bare potential (33), the true vacuum of the system 
is stable. The existence of a negative eigenvalue in actually stable systems is known 
as the problem of a fake instability. To overcome this problem, the valley method 
has been developed recently. B'ElP It is a generalization of the instanton method 
that is based on the valley structure in the configuration space. 



As shown in Figure. 15, an asymmetric bare potential leads to an asymmetric 
effective potential. We show in Figure. 16 results for the energy gap in the cases of 
three values of ho, from bottom to top, ho = 0.02, 0.2, 0.4. For any value of ho, in 
the Ao — > limit, AE approaches y/2. This is because in this limit the asymmetric 
double well approaches a single well. We employ the operator expansion and give the 
truncation N = 12, 14, 16 results. We also plot the results obtained from the valley 
method with fourth and sixth order perturbations. 153 A complementary relation 
between the NPRG and the valley method is observed, just as in the case of the 
symmetric potential. 

As mentioned above, since in the Ao — > limit the potential approaches a single 
well, if we carry out the operator expansion at the potential minimum x = x m i n , the 
NPRG equations never become singular even in the Ao — > region, and we obtain 
m c g ~ y/2. We use this technique for analysis of SUSY QM in the next section. 

We display results for other quantities in Figure. and Fi gurejlq for three values 
of ho, from top to bottom, ho = 0.02, 0.2, 0.4. The expectation value of x, M\ = 
(f2\ x\Q), is shown in Figure. 17, and the variance of x is shown in Figure jl^. The 
NPRG results appear to be perfect on the strong coupling side, while they are 
incorrect in the weak coupling region. 



Change of Potential 

LPAW-H N=16,i„.°- 1 >fy>-0.1 



Energy Gap 
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valley method + 
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valley method + 


6th p.t. 


LPA W-H N=12 




o LPA W-H N=14 




a LPA W-H N=16 




exact result 





Fig. 15. Flow of the potential. 



Fig. 16. Energy gap estimates. 
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Fig. 18. {f2\x 2 \f2) c 



§5. Applications to various quantum systems 

We have seen that the NPRG method is very effective in analyses of non- 
perturbative dynamics in quantum mechanical systems. Here, we apply the NPRG 
method to more non-trivial quantum systems. 

5.1. Super symmetric quantum mechanics 

Here we analyze supersymmetric theory, in which the non-perturbative dynamics 
of the system are crucial. We consider the SUSY QM thepry t _which was introduced 
by Witten as a toy model for dynamical SUSY breaking. Ej)'ES The Hamiltonian is 
given by 



H 



P 2 + W\x) + a 3 



dW(x) 
dx 



1 £>2 

2 



P 2 + V + (x) 

2 



ip 2 + U_(x) 



' tjt2, ^ , 1 dW{x) 



where W(x) is called the SUSY potential. We define the super charges 

Q 1 = ^(a 1 P + a 2 W(x)), 
Q 2 = l -{a 2 P-a 1 W{x)), 

and the Hamiltonian is written as 

H = {4,4} = {Q 2 ,Q 2 }- 

This ensures that the vacuum energy is always non-negative: 



£ = (f2\H\n) = 2 Qt\n) 



> 0. 



(39) 
(40) 

(41) 
(42) 

(43) 

(44) 
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The vacuum energy Eq is the order parameter of dynamical SUSY breaking; 
that is, 

E = => Qi|^) = 0, Q 2 \0) = unbroken SUSY, 
E > => ^ 0, Q 2 \ty ^ broken SUSY. 

Furthermore, the perturbative corrections to are vanishing for any order of the 
perturbation. This is known as the non-renormalization theorem. In fact, with the 
SUSY potential W(x) = gx 2 — x, the potential V+{x) becomes 



V+(x) 



1 o 4 3 f 2 1 
— OX — OX H X + OX . 

2 2 2 



(45) 



The perturbative corrections to the energy spectrum are calculated as 

3 3 
E n = n + -g 2 [2n 2 + 2n + 1] - -# 2 [10n 2 + 2n + 1] 



32 



# 4 [34n 3 + 51n 2 + 59n + 21] + 



(46) 



These corrections to Eq are canceled out at each order of g, and thus there are no 
perturbative corrections. Hence, a non-vanishing Eq is realized only through non- 
perturbative effects caused by the essential singularity at the origin of the coupling 
constant. 



Shape of Bare Potential 
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g-0.1 
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Fig. 19. Bare potentials. 



Fig. 20. Potential flow. 



We analyze thissystem using the LPA W-H equation with iV = 16 for the 
operator expansion. We calculate the effective potential for a wide range of values 
of the parameter g. The case of vanishing g corresponds to the harmonic oscillator 
with a constant term — i, and SUSY does not break in this case (£^o = \ — \ = 0). 
However, SUSY is dynamically broken for any non-vanishing g. Note that for small 
g, the bare potential is an asymmetric double-well, while for g > \/ rhs — 0.31, it 



is 



a single- well, and quantum tunneling is irrelevant (Figure. 19). For any value of g 
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the minimum of the bare potential V+ is at x = 0. Figure. 20 displays the result for 
5=0.24, where the effective potential evolves into a convex form, and its minimum 
turns out to be positive; that is, our NPRG method gives a positive Eq correctly, 
and describes the dynamical SUSY breaking. 



Energy Spectrum 

W(x)=gx z -x , V,(x) 




Vacuum Energy 

W(x).gx a -x , V.(x) 



- exact result 
• LPA W-H N=16 
valley method 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

g 



Fig. 21. Energy spectrum. 



Fig. 22. NPRG and valley method. 



As is shown in Figure|2l], the NPRG results are excellent in the weak coupling 
region and strong coupling region, but not in the region where the bare double- 
well potential becomes deep. In this intermediate region (0.1 < g < 0.2), we cannot 
obtain reliable results because of large numerical errors, while the valley method 
works very well, as shown in Figurej2^. The valley method evaluates the ground 



state energy as Eq = and reproduces the exact value in the weak coupling 

region. 0) However, it does not work in the strong coupling region (g > 0.31), where 
the valley instanton is no longer a good approximate solution of the valley equation. 
Again, we find that the two methods are complementary. 

5.2. Two particle systems 

Next, we apply the NPRG method to quantum many particle systems. As 
the simplest system, we analyze two particle (^1,^2) dynamics with the following 
potential Vo(</?i, V2): 



1 



1 



(47) 



Without the interaction F(<pi,(p<2) between the two particles, the four degenerate 
ground states are mixed by tunneling, splitting into three (fx «-> (p% symmetric states 
and one anti-symmetric state. For the interaction F(ipi, ^2), we now choose ipi *-* ipi 
symmetric interactions and investigate how this interaction affects the energy levels 
of three symmetric states. 
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The LPA W-H equation for two particles is written 

4^4^>4' 1+ ^)' (48> 

where "Tr" represents the trace over the subscripts a, b which correspond to the two 
particles. We consider three types of interactions, 

F((fi, <p 2 ) = C(pnp 2 , C 2 Oi - (p 2 ) 2 , C±(ipx - ip 2 ) 4 - (49) 

In the cases of the second and third types, for C2, C4 > the interaction is attractive, 
and for C 2 , C4 < it is repulsive. Since we now treat only the {ip\ <-> 992) state, it is 
convenient to convert the variables from ((fi,(p 2 ) to (£1,22) as follows: 



(x x \ = wi-i\ (<pi 

\x 2 J 



(50) 



The LPA W-H equation for (xi,X2) has the same form as (|48|). The lowest energy 
splitting for symmetric state, AE = E\ — Eq, is expressed in terms of the effective 
mass of x 2 . Of course, in the C = case, this is equal to the effective mass in one 
particle system. 

The bare potentials are written, corresponding to (149]), as 



V (x) = - C)x\ + ^x\ + + C)x\ + ^x\ + 3A xfx|, (51) 
V (x) = i(-l + 4C 2 )xj + ^xf - \x\ + ^x\ + 3Aox?xi, (52) 
^o(x) = -^1 + (y + 4C 4 ) xl - ^ + ^x\ + 3X xlxl (53) 

We analyze these systems for small C, C 2 and C4. We set Ao = 0.2, which is in the 
parameter region where the NPRG works perfectly in previous analyses. The LPA 
W-H equation was solved numerically using the operator expansion with iV = 12. 
We also calculated AE from the first order perturbation theory with one particle 
Schrodinger wave functions. 



The results for small C,C 2 and C4 are shown in Figures 23, 24 and 25. We 
see that the NPRG results and the Schrodinger wave function results are almost 
the same in these small interaction regions. These results indicate that an attractive 
interaction (C 2 , C4 > 0) causes AE to decrease, and a repulsive interaction (C 2 , C4 < 
0) causes it to increase. 

Here we have shown that for multi-particle systems with interactions, the NPRG 
method can be applied equally without any change of formulation. We are now 
carrying out calculations to obtain non-trivial relations between particle interactions 
and tunneling enhancement /suppression. These results will be reported elsewhere. 

§6. Summary and Outlook 



We have applied the NPRG method to various quantum systems and used it to 
analyze non-perturbative physics. Even in the first stage of approximation, LPA, we 
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Energy Gap 

double-wells (\=0.2) with Cij),^ interaction 



W 0.605 




Fig. 23. Linear interaction. 
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double-wells (X =0.2) with C(0 r ^ 2 ) interaction 



double-wells (X =0.2) with C(0 r 4g interaction 
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»LPA W-H N-12 
- SCH , O(C) 



Fig. 24. Quadratic interaction. 



Fig. 25. Quartic interaction. 



successfully evaluated the non-perturbative quantities that should be given by the 
summation of all orders of the diverging perturbative series. We also found that for 
non-perturbative quantities characterized by an essential singularity, the LPA W-H 
equation again works very well in the region where the instanton-type method breaks 
down, i.e. the strong coupling region. However, NPRG is not effective in the weak 
coupling region, due to large numerical errors. In these regions, the approximation 
used to solve the NPRG equation should be improved in order to obtain correct 
results. To summarize, the NPRG method and the instanton (or valley) method 
play complementary roles. Also, from a practical point of view, the NPRG method 
is a useful new tool for analysis of various quantum systems in a wide parameter 
region. We have obtained good non-perturbative results for SUSY QM. We also 
showed that interacting quantum particles can be treated in a similar way. 

In the flow diagrams, we observed singular behavior in the small coupling region, 
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and found that it becomes more singular under low-order truncation of the operator 
expansion. The origin of the difficulty which we encounter in our NPRG analysis 
resides in the approximation scheme we employed. We must develop 'better' approx- 
imations, which may depend on the individual systems under study. We also need 
to study in detail how to extract physical information from the effective potential 
and the effective action. 



Appendix A 

First Pole Dominance 



Here we confirm the first pole dominance in the two-point function of anharmonic 
oscillators. In the local potential approximation, the two-point function is given by 
following 



E 



Dr. 



LPA 



E 2 + ml 



<•[[■ 



This substitutes one pole for an infinite number of poles. Therefore, if the multi- 
pole contribution becomes significant, the correspondence (24) must be wrong. We 
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Fig. 26. Di for single well potential. 



Fig. 27. D\ for double well potential. 



evaluated the first pole coefficient D\ for a single-well (31) and a double-well (33) by 
solving the Schrodinger equation numerically. The results are displayed in Figures 
26 and ^jj. 

We should note that the relation J2n D n = 1 always holds. For the single- well 
potential, the first pole dominates almost completely. This corresponds to the fact 
that the results obtained with the LPA W-H equation reproduce the correct results. 
On the other hand, for the double-well potential, the first pole dominance begins to 
disappear in the region near Ao = 0.1 — 0.15, where the results obtained with the 
LPA W-H equation become poor. 
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